Complexity of Quantum States and Reversibility of Quantum Motion 
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We present a quantitative analysis of the reversibility properties of classically chaotic quantum 
motion. We analyze the connection between reversibility and the rate at which a quantum state 
acquires a more and more complicated structure in its time evolution. This complexity is character- 
ized by the number M{t) of harmonics of the (initially isotropic, i.e. M{0) = 0) Wigner function, 
which are generated during quantum evolution for the time t. We show that, in contrast to the 
classical exponential increase, this number can grow not faster than linearly and then relate this 
fact with the degree of reversibility of the quantum motion. To explore the reversibility we reverse 
the quantum evolution at some moment T immediately after applying at this moment an instant 
perturbation governed by a strength parameter ^. It follows that there exists a critical perturbation 
strength, ^c ~ y/2/M.{T), below which the initial state is well recovered, whereas reversibility dis- 
appears when ^ > £,c{T). In the classical limit the number of harmonics proliferates exponentially 
with time and the motion becomes practically irreversible. The above results are illustrated in the 
example of the kicked quartic oscillator model. 

PACS numbers: 05.45.Mt, 03.65.Sq, 05.45.Pq 



I. INTRODUCTION 

Strong numerical evidence has been obtained that the 
quantum evolution is very stable, in sharp contrast to 
the extreme sensitivity to initial conditions and rapid loss 
of memory which is the very essence of classical chaos. 
In computer simulations the latter effect leads to practi- 
cal irreversibility of classically chaotic dynamics. Indeed, 
even though the exact equations of motion are reversible, 
any, however small, imprecision such as computer round- 
off errors, is magnified by the exponential instability of 
trajectories to the extent that any memory of the initial 
conditions is effaced and reversibility is destroyed. In 
contrast, almost exact reversion is observed in numerical 
simulations of the quantum motion of classically chaotic 
systems, even in the regime in which statisticalphenom- 
ena such as deterministic diffusion take place [l|]. 

It is intuitive that the physical reasons of this strik- 
ing difference between quantum and classical motion arc 
rooted in the quantization of the phase space in quan- 
tum mechanics. If we consider classical chaotic evolution 
(governed by the Liouville equation) of some phase space 
distribution, smaller and smaller scales are explored ex- 
ponentially fast with time. These fine details of the den- 
sity distribution are lost due to finite accuracy (inevitable 
coarse-graining) in numerical simulations, and therefore 
the reversal of time evolution cannot be carried out. On 
the other hand, in quantum mechanics one expects that 
the structure of a quantal phase-space distribution, e.g. 
of the Wigner function, has resolution limited by the 
size of the Planck's cell. While the mean number of 
Fourier components of the classical phase-space distri- 
bution grows exponentially in time for chaotic motion. 



the number of the components of the Wigner function at 
any given time is related to the degree of excitation of 
the system (see for example eq. (|69|) below) and there- 
fore unrestricted exponential growth of this number is 
not physical Q, y, 3 ■ This fact implies substantially sim- 
pler phase space structure in the case of quantum mo- 
tion as compared with that of classical chaotic dynam- 
ics. We demonstrate below that the mean number of 
Fourier Harmonics is a simple relevant measure of struc- 
tural complexity which in turn is related to fundamental 
properties such as dccohcrencc and entanglement. 

In spite of the above arguments, a rigorous link be- 
tween the intuitively expected different degree of re- 
versibility of quantum and classical motion and the struc- 
ture developed by the phase-space distributions during 
dynamical evolution has never been established. The 
purpose of the present paper is to clarify this problem. 

Following the approach developed in [5| we consider 
first the forward evolution 



m = u{t)muHt) 



(1) 



of an initial (generally mixed) state p(0) up to some time 
t — T. A perturbation P(^) is then applied at this 
time, with the perturbation strength ^. For our pur- 
poses, it will be sufficient to consider unitary perturba- 
tions P(^) = e~*^^, where y is a Hermitian operator. 
The perturbed state 



p{T.o^mf>{T)p\o 



(2) 



is then evolved backward for the time T , thus obtaining 



the reversed state 



mT,0 = UHT)^{T,OU{T) 



(3) 



where P(e,T) = e'^^^f^), with V{T) = U^(T)VU{T) 
being the Heisenberg evolution of the perturbation dur- 
ing the time T. 

Finally, we investigate the distance between the re- 
versed /5(0|T, ^) and the initial /5(0) state, as measured 
by the Peres fidelity p] 



F{C;T) 



Tr[^(0|r,Op(0)] 



TiMt,om] 



Tr[p2(t)] 



Tr[p2(0)] 



^(?;t)lt=T 



(4) 



t=T 



This quantity is bounded in the interval [0, 1] and the 
distance between the initial and the time-reversed state 
is small when F{^]T) is close to one Q- In particular, 
-^(C; T) = 1 when the two states coincide. The second 
line in Eq. Q is a consequence of the unitary time evolu- 
tion and will allow us to relate the distance between the 
initial and the reversed state to the complexity of the 
state p{t) at the reversal time t ~ T. This relation will 
play the key role in our further analysis. 

As we have already mentioned above, we characterize 
the complexity of a quantum state by the structure of its 
Wigner function. The basic idea here is that a quantum 
state is complex if this function has a rich phase space 
structure, which can be naturally measured by the num- 
ber A4{t) of its Fourier harmonics. We will show that 
the fidelity F{^;T) is a decreasing function of the com- 
plexity A4{T) of the state p{T). We will prove further 
that, after the Ehrenfest time scale, namely after a time 
logarithmically short in the effective Planck constant of 
the system, the number A4 (t) of harmonics increases not- 
faster than linearly with time. On the other hand, in clas- 
sical chaotic dynamics the number of harmonics A4c{t) 
of the classical phase-space distribution function grows 
exponentially in time. We will then ascertain that the 
initial state is well recovered as long as the perturbation 
strength ^ is much smaller than a critical perturbation 
strength ^^(T) - 1/M{T). Therefore, Cc(T) drops expo- 
nentially with T in the classical case and not faster than 
linearly for quantum evolution. This fact explains the 
much weaker sensitivity of quantum dynamics to pertur- 
bations as compared to the classical chaotic motion. 

The paper is organized as follows. In Sec. II, we review 
the main concepts and definitions of quantum mechan- 
ics in phase space relevant for our work. In Sec. Ill, we 
discuss the kicked quartic oscillator model, used in the re- 
maining part of the paper as a test bed to illustrate the 
stability properties of classically chaotic quantum mo- 
tion. In Sec IV, the evolution in time of the harmonics 
of the Wigner function is studied in detail. In particular, 
the complexity of the quantum state p{t) is quantified by 



looking at the sensitivity of the system to an infinitesi- 
mal perturbation. In Sec. V, the degree of reversibility of 
quantum motion, as measured by the Peres fidelity, is re- 
lated to the number of harmonics developed by dynamics 
at the reversal time t = T. A more detailed study of the 
reversibility properties of motion is carried out in Sec. VI, 
where we investigate the properties of the time-reversed 
state by studying the harmonics of its Wigner function. 
Finally, the main results of our paper are summarized in 
Sec. VI. 



II. QUANTUM DYNAMICS IN THE PHASE 
SPACE 

The phase-space representation of quantum mechan- 
ics is a very enlightening approach which allows a direct 
comparison with classical mechanics 0, 0] . In this sec- 
tion, we briefly review the main aspects of the phase- 
space approach which are relevant for our work. 



A. The Wigner Function 

Let us consider a nonlinear system whose dynamics is 
governed by the Hamiltonian operator H = H{a\ a; t) 



HW^f 



7t, 



;) + ij(^' (a^ , a; t) with the time-independent 



unperturbed part H^^^^ which has a discrete energy spec- 
trum bounded from below so that we can assume that 
all its eigenvalues En ^ 0. Here a^,a are the bosonic, 
[a, d^] = 1 , creation-annihilation operators. We will use 
the method of c-number a-phase space borrowed from 
the quantum optics (see for example (lOl . [III). This 
method is equally suitable for analyzing both the quan- 
tum and classical evolutions. It is, basically, built upon 
the basis of the coherent states |a). The latter are de- 
fined by the eigenvalue problem a\a) = -j=\a), where 
a is a complex variable independent of h. An arbi- 
trary coherent state is obtained from the ground state, 

\a) ~ D (-j=] \0), with the help of the unitary displace- 
ment operator D (A) = cxp(A d^^ — A*d). 

Using the displacement operator D as a kernel, one can 
represent any operator function G{t) = G(d^j_a; t) in the 
form of the operator Fourier transformation [ll| 

Gia\a;t) ^ ^ J d^i^Giv* ,fj;t) D{ri) , (5) 

where G{r]* , 77; t) is a numerical function of two indepen- 
dent complex variables 77* , 77 and the integration runs over 
the complex 77-planc. The inverse transformation 



G(r/*, 77; i) =Tr G{a\a-t) D\ri) 



it^ 



(6) 



is immediately obtained with the help of the orthogonal- 
ity condition for the displacement operators. 



1 



■Tr 



b\ii)b{'n) =,5(2)(r/-r;). 



(7) 



By using transformation ([5]), the standard quantum- 
mechanical formula {Q) = Tr p(t) Q for the mean ex- 
pectation value (Q) of a dynamical variable Q [Q is rep- 
resented by the operator Q = Q{a^,a)] in a generally 
mixed state p{t) = p{a\a\t) can be written in a way 
formally equivalent to the classical phase space average: 



J (faW{a\a;t)Q{a*,a) . 



(8) 



The final form is readily obtained after defining the 
Wigner function W{a*,a;t) and Q{a*,a) as c-number 
Fourier transformations of p{rf ^rj\t) and Q(r]*,r]), re- 
spectively. The Wigner function in the a-phase plane 
is connected to the density operator p{t) as 

Wia*,a;t) ^ ^J £rj cxp{jj^^ r^*^) p{rj*,r^;t) ^ 

^Jd'rj exp (r;* f^-V^) Tr [p(t) D{r^) 

(9) 
(this corresponds to the Weyl's ordering of the creation- 
annihilation operators). Similarly, 

Q{a*,a) = i/ A cxp (rj^-rj*f^) Qir^^-q) = 
i/Aexp(7r^-77^)Tr[Q^(7,)" . 

(10) 
It follows from the definition ^ that the Wigner func- 
tion is normalized to unity: 



/ 



d?aW{a\a]t) =Trp{t) = 1. 



(11) 



The Wigner function is real but, unlike its classical coun- 
terpart, is not in general positive definite. 

With the help of the Wigner function the Peres fidelity 
(|3]) can be expressed as 



F{t,T) 



j ffiaW {a\a]Q)W {a\a: 



o|r,c) 



^ (PaW {a\a]T)W {a* ,a;T,C) 



(12) 



The important advantage of this representation is that 
it remains valid in the classical case when the Wigner 
function reduces to the classical distribution function, 

M/(=)(a*,a;i). 



where a = ^/l e~^^ . Then the normalization condition 
simply implies that J„ (i/Wo(/;i) = 1. There arc no 
restrictions on W™ when m ^ Q. 

The amplitudes Wm can be expressed in terms of the 
matrix elements {n + m\p\n) along TTi-th subdiagonal of 
the density matrix. Indeed, using the well-known |12l | 
matrix elements of the displacement operator in the basis 
of the eigenvectors \n) of the unperturbed Hamiltonian 



{n + m\D{vi)\n) 



[n + to) 



ry-V-^I^I^L^Td^l'), (14) 



(n, TO > 0) where L'^{x) is the Laguerre polynomial, the 
77-integration in the second line of Eq. ^ can be carried 
out explicitly. We obtain finally 



M^™(/;i) = te-*^Er=o(-irv(^^T;^- 

{AI/h)^Ll^{AI/h){n + m\p{t)W) 



(15) 



w:„>o- 



when m > and W- 

With the help of the orthogonality condition for the 
Laguerre polynomials Eq. (|15p can be inverted, thus ob- 
taining 



(n + TO|p(i)|n) = (-l)"2yC|Erx 



(16) 



C. Time - Evolution of the Wigner Function 



The Wigner function satisfies the evolution equation 

i^^Wia*,a-t)=CgWia*,a;t). (17) 

Here Cqis the Hermitian "quantum Liouville operator" 
Cq whose explicit form is obtained by mapping the stan- 
uation 



dard equation 



4,m = \i^^m] 



(18) 



onto the a-phase space. Supposing that the Hamiltonian 
can be presented as a Hermitian sum (finite or infinite) of 
products of the creation-annihilation operators and using 
then the definitions of the phase space images, Eqs. ([9]), 
([TU)) . it is possible to show that (see Ref. |11| ) 



B. Harmonics of the Wigner EYinction 

We define the harmonic's amplitudes Wm{I',t) of the 
Wigner function by the Fourier expansion 



_. 00 
W{a*,a;t)^- V W^{I;t), 

77- / J 



(13) 



Ca 



Ha 



^^,a+'^-^ 



2 Sq' 



2 da' 



rt [a + 9 Q„ , ck „ Q » 



2 da' 



2 da' 



(19) 



where Ti.{a* , a) is the phase-space image of the system's 
Hamiltonian. The arrows above the derivatives mean 
that they act only on the arguments of the Wigner func- 
tion but ignore the (a*, a)-dependcncc of the phase-space 



operator TC itself. In the classical limit h = the operator 



Cq has the standard classical form 



Cc 



dHc{a* ,a;t) d 
da da* 



dHc{a* ,a;t) d 
da* 9a' 



(20) 



where the classical Hamiltonian function coincides 
with the diagonal matrix element, Hc{a* ,a;t) = 
{a\H'^^^a^ , a)\a) of the normal form H^^^ of the quan- 
tum Hamiltonian operator. In other words, this function 
is obtained from the quantum Hamiltonian by substitut- 
ing a -^ a/V", o,^ — >■ a* I \lh . 

The outlined phase-space app roach is quite general and 
can be readily extended [10, lll| to systems with arbitrary 
number of degrees of freedom and is applicable to any sys- 
tem whose Hamiltonian can be expressed in terms of a 
set of the bosonic creation-annihilation operators [I3 • In 
this approach, the classical distribution function as well 
as its quantum counterpart are described in terms of the 
same phase-space variables, thus allowing a straightfor- 
ward comparison of the two dynamics [3|, [ 



III. THE MODEL 

In order to discuss the reversibility/complexity prop- 
erties of quantum motion, we consider, as an illustrative 
example, the kicked quartic oscillator model, described 
by the Hamiltonian [3, [ij, llJ] 

ii = H(a^ a) = hujQh + h^n^ - Vhg{t){a + a^), (21) 



a) a, 



1. 



where g{t) ~ goj^s^i^ ~ ^)' ^^^'^ '"^ 
In our units, the time and parameters /i, cjq as well as 
the strength of the driving force are dimensionless. The 
period of the driving force g(t) is set to one. The cor- 
responding classical Hamiltonian function He can be ex- 
pressed in terms of complex canonical variables a, ia* 
which are related to the classical action-angle variables 
/, 9 via a = ^/le~^^ ,a* — ^/le^^ . It reads 



He = ujo\a[ 



■ g{t){a* +a). 



(22) 



Detailed analytical semiclassical analysis of the quantum 
motion of the model (PT|) has been presented in [ij, \lM ■ 
The quantum Liouville operator for the model (I2ip can 
be derived from Eq. (J19p and reads 






UJQ 



-h 

r*{kick) 



2 da* da 



~9it) 



f*{kick) 

-2|ap 

d 
da* 



a- 



d_ 
da 



(23) 

The essential difference between this operator and the 
corresponding classical Liouville operator Cc is the pres- 
ence in Eq. (P5|) of a term proportional to T? which con- 
tains a second derivative over the phase space variables. 



This term drastically changes the spectrum of the unper- 
turbed part Cq of the Liouville operator and thereby 
modifies the evolution of the Wigner distribution func- 
tion W with respect to the classical distribution function 
W^. Indeed, while in the classical limit h = Q the factor 
tjQ + 2|ap ~ ijJo + 2I = ujc{I) is the continuous frequency 
of the classical quartic oscillator, the quantum operator 
K = —■kfi^-oAra — l-2|ap has a discrete spectrum: consid- 



ering the real and imaginary parts {ai, a2) of the variable 
a as cartesian coordinates, we obtain 



K 



d^ 
da\ 



da\ 



^2{a\ + a\). (24) 



Therefore, the operator K is formally equivalent to the 
Hamiltonian operator of a two-dimensional isotropic os- 
cillator, with the frequency v = \ and the mass ^ = 4. 
After introducing the standard annihilation-creation op- 
erators 



Al,2 

^1,2 



ai,2 



ai,2 



1 d 

VsK 9q!i,2 ' 



(25) 



this operator transforms into K = HN, where N = Ni + 
N2 ~ ^1^1 +A2A2 is the operator representing the total 
number of fictitious quanta linearly polarized in the a- 
plane. The operator 



M 



, d 
da* 



d_ 
da 



.d_ 

'de 



(26) 



is proportional to the angular momentum projection op- 
erator along the axis orthogonal to the a-plane and has 
the discrete eigenvalue spectrum m ~ 0, ±1,±2, .... We 
can therefore conclude that the spectrum of the operator 



^(0) 



(uQ + hN) 



M 



(27) 



is also discrete with eigenvalues (Ao)n.,„ ~ {ujq + hn)m. 
Since the operators N and M commute, operator ([?7|) is 
Hermitian. 

The two linearly polarized quanta introduced above 
are coupled to each other because the angular momen- 
tum operator hM = i I A]^A2 — AjAi I is not diagonal 

in the chosen representation. However both the opera- 
tors hN and hM are simultaneously diagonalized after 
introducing the operators 



A, 



1 

71' 



^/h d 
Tda* 



A_ 



1 



Vh d 

'Yda 



(28) 



which describe the circularly polarized quanta. In the 
new representation 



N = N+ + N^, M = {N+ - 7V_), 



(29) 



where N± ~ A'j. A± are the operators representing the 
numbers n± of circularly polarized quanta. These new 



quanta arc decoupled: 

4°' = (^0 + ^iv) M = (£o)+ - (Co). 



(30) 



The eigenvalues (Ao ),!+,«_ 



"+ ^ " of the operator 
([50)1 are determined by the distances between the un- 
perturbed energy levels En^ — Huiq n± + h^ n\ corre- 
sponding to the excitation numbers n± = 0, 1,2, .... The 
representation (|30|1 is the most convenient for numerical 
simulations. 

The driving perturbation also decouples in the chosen 
representation: 



£{kick) _ Mkick) _ Mkick) 



(31) 



It then immediately follows that the one-period unitary 
evolution operator J- for the Wigner function gets fac- 
torized as 



T = T+P_, T±= e-*(^°)± b± te) 
^±fe)=cxp(z^(i±+in). 



0^( 



(32) 



The complete set of the eigenvectors |n_(.n_) = 

|n+)|n_) of the unperturbed operator tq constitutes 
the excitation number reference basis for the density ma- 
trix, 

p{n+,n_;t) = (n+|/5(i)|n_) . 
The evolution from time t to time t -I- 1 reads 
p(n+,?i_;i+l) = {n+\Uip{t)Ul\n^) = 



Y.n'n' {n+\T+\n\)p{n'+,n'_-t){n'_\TAn^Y 



Here 



^7l = [/(i=l) = e-i^""^(^^ 



(34) 



(35) 



matrix of a finite size A''. The main disadvantage of this 
approach is the violation of unitarity of the Floquet op- 
erator: its norm is not conserved when the size of a quan- 
tum state becomes of the order or larger than N . 

The second approach is based on the truncation of the 
Hermitian matrix X„„' — (fi|-% (a + a^) |?^')- Then we 

numerically diagonalize this matrix, X = vXdv\ with 
the aid of a unitary matrix v. In such a way we obtain 
the truncated kick matrix D{i-yp) = vexp{iXd)v^ in the 
excitation number basis. Multiplying it by the, diagonal 
in this basis, free rotation matrix exp ( — ^7?*^°-* j (trun- 
cated to the same size N) and finally diagonalizing the 
obtained unitary matrix, we arrive at a truncated ap- 
proximation F ~ VFiiV^ of the Floquet operator ([551) 
in the excitation number representation. The price paid 
is the artificial boundary condition at n = A^ which in- 
fluences the evolution when the size of a quantum state 
becomes close to that of the truncated Floquet operator. 
Even though the final diagonalization problem is, by it- 
self, rather time-consuming the great advantage of this 
approach is that the computation time does not depend 
on the considered^ duration of the^ evolution. This is due 
to the fact that U{t) = F* = VFl^V^ and therefore, in- 
dependentlj of t, it is sufficient to multiply 3 matrices to 
construct U{t). 



IV. HARMONICS DYNAMICS 

In this section, we study the time evolution of the har- 
monics of the Wigner function for the kicked quartic os- 
cillator model. The Wigner function of the initial state is 
taken to be isotropic, so that only the zero angular har- 
monic Wq of the Wigner function is different from zero. 
In the course of dynamical evolution the Wigner func- 
tion becomes more and more structured and this reflects 
in the excitation of a growing number of harmonics. The 
complexity of the Wigner function at time t is measured 
by its sensitivity to an infinitesimal perturbation. 



Initial Condition 



is the (Floquet) operator for our model, namely the one- 
period unitary evolution induced by the Hamiltonian 

(ED). 

Our numerical simulations are based on the combined 
application of Eqs. ([4]), ([9]), and (p4|) . We calculate nu- 
merically the truncated Floquet matrix Fnn' in the ex- 
citation number representation. Two different strategies 
have been used. 

The first approach is based on Eq. (fT4|) . which relates 



matrix elements of the kick operator D ii-^\ to the La- 

guerre polynomials. The free rotation operator 6""^^ is 
diagonal and its calculation is trivial. Then the product 
of the matrices e~i^ and D is truncated to a square 



We choose the initial state to be an isotropic mixture 
of coherent states: 

/oo 
d2a7'(|a|2)|a)(a| = ^p„|7i)(n|, (36) 
n 



where 



P- = -^j^ dIViI)e-'/^[l/hj , I^\a\\ (37) 

Here and in the following, a circle above a dynamical 
variable denotes its value at the time t — 0. To simplify 
further analytical considerations we suppose a Poissonian 



100 




FIG. 1: Mean value ({/)t — ^)/go as a function of time t. 
Squares and triangles correspond to (ft, (/o) = (l, 2) and (2,3), 
full and empty symbols to A = and 50. The straight dashed 
line coresponds to the classical diffusion law {I)t = A + g^t. 



.7;^ ] in the exci- 
tation number space, which implies the exponential form 

O O 

V{I) — -:^ e~^' ^ of the distribution in the space of co- 
herent states and, correspondingly, the isotropic Gaus- 
sian initial Wigner function 



W{a\a]0) 



1 



li + h/2 



(38) 



In the particular case of the pure ground state, A = 
0, the Wigner function occupies the minimal quantum 
cell with the area h/2 . It is worth noting that this area 
would be twice as much in the case of the normal ordering 
of the creation annihilation operators (Husimi function) 

©[HI. 

For the initial conditions (|36p the classical dynamics of 
the model (|2ip becomes chaotic when the kick strength 
parameter go exceeds a critical value go.c ~ 1- The an- 
gular phase correlations decay exponentially (we have 
checked this fact numerically) and the mean action grows 
diffusively with the diffusion coefficient D = g^. Our nu- 
merical data presented in Fig. [1] demonstrate the corre- 
sponding "quantum diffusion" phenomenon described in 
the following subsection: (/)t = A -I- gg t as in the clas- 
sical case, until a time t* after which the quantum to 
classical correspondence breaks down 0, [ij] . 



B. Zero- Harmonic Evolution and Quantum 
Diffusion 



According to Eq. ((T0|), the zero amplitude Wo(/;i) is 
determined as 



Wo[I;t)^-e- 



\'Y.i-irU,{U/K){n\p{t)\n) (39) 



n=0 



by the diagonal matrix elements of the density operator. 
Our numerical simulations (see Fig. ^ show that these 



matrix elements decay along the main diagonal on aver- 
age exponentially at any moment t after a short initial 
interval (see below). Neglecting fluctuations we therefore 
assume that 



Pnnit) = ("-|/0(i)l"-) 



2 



[n\U{t)\a) 



Tn^^oPn' {n\U{t)\n') =[l-e 



Mm 



-y{t)n 



(40) 

with a function v{t) which depends on the system's dy- 
namics. The approximation is similar to coarse graining 
of the classical distribution function. Correspondingly, 



Wo{I;t) 



2 , fv 

— tanh — 1 cxp 

n \2' 



2 , /i^ 

— tanh — 

h \2 



I 



(41) 



The dependence on time of the mean value of a func- 
tion Q{I) of the action operator I = tin is then computed 
with the help of the formula 



(g>t= / diQ{i)Womt) 

Jo 



(42) 



where the phase-space image of the operator Q is calcu- 
lated similarly to the zero harmonic amplitude Wo{I;t) 
of the Wigner function (see Eq. (fTO|l ). For example the 
image of the operator G{k) — e~'^^ , which generates the 
moments of the amplitude WQ{I;t) is easily found to be 



G{k-I) 



oKS/2 



cos\i{tih/2) 



exp 



-2tanh(Kft/2)- 

h 



(43) 



The mean value of the generating function of the action 
momenta is readily obtained from eq. (|42p 



{G{K-I)), = e' 



sinh(i^/2) 



sinh(://2 4- Kh/2) 



(44) 



In particular, the mean action equals {I)t = 
-d{G{K-I))t/dK\^^^ = I [coth(f) - 1]. This formula 
relates the function v{t) to the evolution of the action, 
g-i^(t) _ (/)j/((/)(-|-;i)^ and allows us to represent finally 
the coarse-grained amplitude of the zero harmonic in the 
form 



»•„(/;() = 



(7)7T|°"= 



/ 



{I)t 



(45) 



The time dependence of the mean action (/) t is shown in 
Fig.[T] 



C. Evolution of Nonzero Momenta: Complexity of 
Quantum States 

The paramount property of classical dynamical chaos 
is the exponentially fast structuring of the system's phase 
space on finer and finer scales. In particular, the num- 
ber JV[[t) of angular harmonics, that is the number 
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FIG. 2: (color online) Distribution of the diagonal elements 
pnn{i) of the density matrix versus n, at h = 1, go ~ 2. 
Left panel (a): mixed state (A = 25) and, from bottom to 
top, t = 10, 30, and 50 (these curves are scaled by factors 
0.01, 0.1 and 1, respectively). Right panel (b): pure state 
(A = 0), i = 10 (bottom, scaled by a factor 0.01) and t = 50 
(top). Straight lines show exponential fits, corresponding to 
the coarse-grained distribution (|40p . 



of appreciably large harmonic's amplitudes VF'-^' m(-^;0 
in the Fourier expansion p^ of the classical distribu- 
tion function W^'^^{a* ,a;t) grows exponentially in time. 
Namely A4{t) ex e*/^=, where the characteristic time Tc 
goes to infinite when the classical Lyapunov exponent 
vanishes. A simple consideration shows that the expo- 
nential regime cannot last long in the case of quantum 
dynamics. Indeed, in the terms of our auxiliary two- 
dimensional linear oscillator where the functions e""'' 
are the eigenstates of the operator M defined in p6|) 
the mean number of harmonics A4{t) ^ {\n^ — n_|)t < 
{N)t = (n+ + n-)t ~ {I)t/h. Therefore the exponen- 
tial upgrowth is possible only for 6*'"^= < {I)t/h, namely 
for a time t < ts ~ Tcln -^-jp . Since the mean action in- 
creases only linearly in time, tE is basically the Ehrenfest 
time |l3| , logarithmically short in h. 

In order to ascertain how complex the quantum state 
became by the time t we use as a probe a phase plane 
rotation by the angle 59 = £,. Such a rotation is gener- 
ated by the unitary transformation P{£,) = e~*^" of the 
density matrix p{t). The effect of such perturbation, as 
mentioned above, is characterized by the Peres fidelity 

O- 

From the second line of Eq. Q we obtain 



Fi^;t)^ J2 cos[e(n'-n)] 



Ami 



n.7i'—0 



EZoikp'im) 



(46) 



The diagonal {n' = n) contribution 



F,{t) = £ 



{n\p{t)\n) 



',T.z,wm) 



(47) 



does not depend on ^. Taking into account that F{^ ~ 
0;i) = 1, we relate this contribution to the non-diagonal 



part of the density operator: 

oo 

Fait) = 1-2 Y. 



{n'\p{t)\n) 



(48) 



After substitution of this expression into (|46p we obtain 



F{^-t) = 1 - 2 ^ sin^ (^m/2) W™(t), 



(49) 



where 



Since 



°° \{n A~ m\p{t)\n) 



i^oT.n'=Mp'm^') 



2Em=iEr=o {n + m\p{t)\n) 
we can convert Eq. ([501) ™to 
>V„.(i) 



(50) 



(51) 



(2-5,„o)i:~=o 



(n+?n p{t) n) 



i+2 2:™=iEr=o 



/e~/=oP' '.'(*) 



(?i+m p{t) \n) 



(52) 



^,7 = P' '„'(*) 



where Pnn{t) is a shorthand notation for (ri/5(i)m). 

Equivalently. one can express W,„(i) in terms of the 
amplitudes p5)) of the Wigner function: 



(53) 



(54) 



W™(i) = 

The completeness condition 

|(4//;i)-"e4^/'"^5(/-/') 



has been taken into account in deriving formula (|53p . 
Similarly to Eq. (|12p . expression (j53p remains valid in 
the classical limit, provided that the harmonics of the 
classical distribution function are used. 

Since the normalization condition X]m^'"(0 = 1 
holds, the quantities Wmit), m > give the probabil- 
ity distribution over the harmonic's numbers m. Now, in 
the spirit of the linear response theory, we consider an 
infinitesimally small rotation angle ^ ^ and hold only 
the linear term of the power expansion of the density 
operator. The Eq. (PO)) reduces then to 



F(e;t)«l 



\e{m^)t. 



?=o 



E™=i"^'W™W 



(55) 



Numerical simulations (see Fig. [3]) show that, apart from 
small fluctuations, the distribution Wm(t) decreases with 
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FIG. 3: Distribution of harmonics VVm(i) as a function of m, 
with parameter values as in Fig. [S] 



m monotonically and exponentially. Therefore the quan- 
tity \J {m?)^ gives an estimate of the number M. of har- 
monics developed up to time t and can be considered as 
a suitable measure of complexity of the Wigncr function 
at the time t. 

According to the first line of Eq. @ . the same fidelity 
can, alternatively, be presented in the form 



ramctcr and investigate the differences between the ini- 
tial and reversed states as a function of ^. To this end we 
first analyze the fidelity F{^;t) at an arbitrary moment 
t as a function of ^. 

Using Eq. (|1^ . we can expand F(^, t) in terms of the 



even moments [r 
tribution ((501) 



■,2fc\ 



(fc = 1, 2, ...) of the probability dis- 



fe=i 



fc+i 



t2k 



(2fc)! 



{m'')t. 



(60) 



We remark that, due to the exponential decay (as a func- 
tion of m) of the distribution Wm{t), all moments {m^^)t 
are finite. 

On the other hand, the expansion of the general ex- 
pression (j56p over the parameter ^ contains on- and off- 
diagonal matrix elements of the powers of the operator 
I{t) = hfi{t) and therefore connect the Peres fidelity to 
the action evolution. Both the equivalent representations 
([50)1 and ([50)) will be exploited below. 



F{^.t) 



EOO PnPn' 

n.n/—() "^^^^ ^ 



i-e \x2it)-2^{t) 



-i(,n{t) I 



Oi^) 



(56) 



where the quantity 



pI 



X2{t) = E ^^^V ((^1"' Wl") - {nHt)\nr) (57) 

is the weighted-mean value of the standard devia tion of 
the excitation numbers at the moment t whereas (j){i) is 
similar wcightcd-mcan value of the cumulative contribu- 
tion 



0(t)=£(l + ft/A)- 



i{t)\, 



(58) 



of the off-diagonal matrix elements. Comparing now the 
^^-tcrms in the both equivalent representations ()55p and 
(|56p of the fidelity we arrive at the following significant 
exact relation between the time behavior of the mean 
number of harmonics (which characterizes the complex- 
ity) on the one hand and of the excitation numbers on 
the other hand 



("^')* 



2 X2(t)-20(t) 



(59) 



The negative second contribution which appears only in 
the case of mixed initial states reduces the number of 
harmonics. 



QUANTUM REVERSIBILITY: THE 
FIDELITY 



To explore the reversibility of the quantum dynamics 
we consider now the perturbation angle ^ as a free pa- 



A. Pure Coherent Initial State 

The theoretical analysis is especially easy to carry out 
in the simple case of the (pure) ground initial state /?„ ~ 
6no- The expression ()56|) reduces then to 



F{tt)^\{o\p{ttm\ = 



i^n{t) 



|0) 



(61) 



This specific initial state is the isotropic coherent state 
|a = 0) which corresponds to the distribution Pdap) ~ 

S^^''{a) in Eq. ([36)) . After first few kicks, a state of 
practically general form is produced. 

Making use of the cumulant expansion we obtain 



F(C; t) = exp 






1=1 



(20! 



(62) 



where the cumulants (connected momenta) are 

X2(t) = (0|(n(t)-(0|n(i)|0))^|0), 

X4{t)^{o\{h{t)-{oHtmf\o)- 

)\{h{t)^{0\n{tmf\0) 



(63) 



and so on. Correspondingly, L{^;t) = —\nF{^;t) = 
^^X2{t) — Y2'?^X4(0 + ■•■ • At a given time t we can retain 
only the lowest cumulant, L(^; t) « S,^X2{t), as long as the 
perturbation strength ^ < 2\/3x2(0/X4(0- However, as 
shown in Fig. [H this approximation fails for larger values 

on- 

To go beyond such a restricted range of values of the 
parameter ^ we observe that the amplitude /(^; t) = 




FIG. 4: (color online) Fidelity F{£^, t) versus perturbation 
strength ^, at h = 0.25, po = 2 at three different times: 
t = 5, Vx2(i) = 80 (squares); i = 25, Vx2(t) = 384 (cir- 
cles); f = 75, y/x2{t) = 1010 (triangles). The dotted and 
dashed curves show decay of the r.h.s. in Eq. (|62|l with, re- 
spectively, only the lowest and the two lowest terms of the 
cumulant expansion being kept. The full curve corresponds 
to the theoretical prediction of the second line in Eq. (|66p . 



(0|. 



,-i^n{t) 



|0) can be readily represented as 



fi^;t) 



n=0 



Wn{t)e' 



i^n 



(64) 



where w„(i) = {n\p{t)\n) = \{n\tj{t)\Q\ is the excitation 
number probabiUty distribution. 

The probabiUty Wn(t) exhibits larger fluctuations as 
a function of n than those in the case of broad initial 
mixtures (compare the two panels of Fig. [5] obtained 
for pure (right) and mixed (left) initial distributions). 
However, at any given time t larger than the Ehrcnfcst 
time Wn{t) decays, on average, exponentially. Assum- 
ing, similarly to (j40p . the exponential ansatz Wn{t) ~ 
[l - e'^(*)] e-^(*)", we obtain 



Wn(t) 



{n)t + 1 



{n)t + 1 



(65) 



where {n)t is the mean excitation number at the moment 
t. The approximation ((65l) leads to the result 



F{^;t) 



l+4(n)t((n)t + l) sin^(C/2) 



1+4x2(4) sin^(?/2) 



(66) 



In the second line we have taken into account that ap- 
proximation (|65p implies the following relation: 

X2it) = {0\n^\0)t - (0|n|0)2 « {n)t {{n)t + 1) . (67) 

(Let us remind in this connection that we have chosen 

a = 0). Opposite to the exact relation X2{t) = 5(™^)t 
(see eq. ([5^) the relation eq. ([^ is vahd only after 
the Ehrcnfcst time i^. Notice that accordingly to the 
first line in Eq. §^ F{2n;t) = F{0;t) = 1, as expected 



since rotation by the angle ^ = 27r is just the identity 
operation. 

Fig. [4] shows that the fidelity decay is nicely described 
by the analytical formula (|66p . It is also clearly seen that 
the fidelity almost vanishes already at very small values of 
the parameter ^, so that the approximation given in the 
second line of Eq. (pS)) works quite well. We can therefore 
conclude that in the considered case of the pure ground 
initial state only the lowest cumulant X2 determines the 
overall decay of the fidelity, whereas the higher ones are 
responsible for the fluctuations in the decay law. 

Using the above exponential ansatz, the expression 
((50)) for the probability distribution Wm{t) simplifies to 



Wm{t) = (2 - S„io) J2'^=0 Wn+m{t) W„ (f ) 



2~S„,o 
2{n)t + l 



{n)t 
{n)t+l 



(68) 



so that this distribution decays with the same slope X(t) 
as the distribution of the excitation numbers (|65|) . To- 
gether with Eq. ([67ll this yields the following relation: 

{m^)t^2x2(t)^2(n)t{{n)t + l). (69) 

Therefore, after an initial interval of order of the Ehrcn- 
fcst time, the number of harmonics of the Wigner func- 
tion increases with time in the same manner as the exci- 
tation number, not faster than linearly. Notice also that, 
as it should be, substitution of the expression (|^ in the 
general formula (|49|) leads again to the same result ((66|) . 
The relatively slow dependence of F{£^; t) and the num- 
ber of harmonics M{t) = ^J {'m?)t on time, which follows 
from expressions ((66)) and ((69)) should be juxtaposed with 
the classical behavior dictated by the exponential insta- 
bility of the classical dynamics. The latter manifests it- 
self, in particular, in the exponential growth of the num- 
ber of harmonics of the classical phase-space distribution 
function Wc{a* ^a^t). To accomplish such a comparison 
we solve the classical Liouville equation with the initial 

phase space distribution Wc{a* , a; 0) = | e ~ of size 5 
which coincides, for a given value of fi,, with the size h/2 
of the Wigner function corresponding to the initial quan- 
tum ground state /5(0) = |0) (0| . The quantum to classical 
transition is explored by keeping 5 constant and consider- 
ing, for smaller and smaller values of h^ initial incoherent 
mixtures ((55)) of size 5 = IS. -\- h/2. A numerical illus- 
tration of such a procedure is presented in Fig. [5| The 
exponential increase of {m?)t up to the Ehrenfest time 
is clearly seen. After that time, a much slower power- 
law increase follows, in accordance, for pure states,with 
the relation ((69)) . Such a behavior is consistent with the 
findings reported in the Refs. @, |3l- 

The first relation in Eq. ()69p allows us to directly con- 
nect the Peres fidelity after the Ehrenfest time with the, 
characterized by the mean number ^J {vi^)t of harmonics, 
complexity of the Wigner function at time t: 
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FIG. 5: (color online) {jin?)t versus t, at go = 1-5, 5 — 0.5. 
Squares, diamonds and triangles correspond to h — 0.01,0.1 
and 1. In this latter case, the initial condition corresponds 
to the ground state, /3(0) — |0){0|. Empty circles refer to 
classical dynamics and the dashed line is an exponential fit to 
these data, {ni?)t ~ exp(a + ht), with a ~ —2.9 and h ~ 4.4. 
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FIG. 6: Distribution Wm of harmonics at times smaller than 
the Ehrenfest time, for parameter values h = 0.01, go = 1.5. 
Data at i = 1 and i = 3 are scaled by a factor 0.1. 



More than that, comparing the terms of the two expan- 
sions, (|60p and ([70]), the former of which is exact for any 
time including times shorter or of the order of the Ehren- 
fest time, we see that also the latter would be always 
correct if 



VI. QUANTUM REVERSIBILITY: THE 
TIME-REVERSED STATE. 

In this section, we study in detail the phase-space 
structure of the time-reversed state characterized by the 
harmonics content of the Wigner function in dependence 
on the perturbation strength ^ and the reversal time T. 

Generally, the Peres fidelity F{^;T), as it appeares in 
the first line of Eq. p2p . is sensitive only to the distortion 
of the zero harmonic of the reversed Wigner function or, 
equivalently, to the redistribution of the excitation num- 
bers. Utilizing the first line of Eq. Q we obtain in the 
terms of the density matrix 



2Er=oe-*" ('^lICoKoh), (A » h), 
where only the diagonal matrix elements 

{n\mT,0\n)^En'=oPn'\{n\P{^,T)\n') 
[n\P{C,T)\n) ' 



(72) 



(73) 



E:^=o(l-'^-')Pn' (n| P(C,T)-1 |n' 



of the reversed density matrix are present. Only the term 
which stays in the second line of this equation remains 
in the limit ^ ^ 0. This limit is shown with the dotted 
line in the Fig. 8. 

The effect of perturbation shows up first in the second 
order with respect to the perturbation parameter ^. Ex- 



panding in Eq. ([751) the operator P{(,, T) 
to the correction of this order we find 

(n|3(0|T,e)|n)« 
p^[l-e{{n\h\T)\n)^{n\h{T)\n)^)] + 

eY.n'^nPn'\{n'\n{T)\n)\^ 



-tin{T) 



up 



(74) 



which leads to Eq. ([56]) again. 

The fidelity ([7^ consists of two different contributions 



r(") 



F„(^;T) -HF("')(^;T) one of which 
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{n\Pi^,T)\n) 



(75) 



{m''\ = ^((m^) 



(71) 



is the weighted mean of pure state fidelities, whereas the 
off-diagonal contribution 



Such a relation is characteristic of the exponentially de- 
caying distribution. Numerical data presented in Fig. [6| 
support such a conjecture. It follows then that in accor- 
dance with the different growth of the number of har- 
monics before and after the Ehrenfest time (see Fig. [5]), 
the slope of the m-dependence of the distribution yV,„ (i) 
is drastically different inside and outside the Ehrenfest 
time scale. This slope decreases exponentially with t in 
the first case and not faster than linearly in the second. 



i^(™)(^;T) 



2Er=o 



I2T=oP', 



■E° 



\H,T)\ 



(76) 
is specific for mixed initial states. 

However, the important information on the harmon- 
ics that survived the backward evolution is absent in the 
fidelity (|72p . To find their number and the correspond- 
ing distribution Wm (0| ^; T) we perturb the reversed den- 
sity matrix /5(0|T, ^) by means of the probing operation 



11 



P{(,') = e~'^ ", with a new infinitcsimally small rotation 
angle ^' . In the same manner as in Sec. IIVCI we obtain 
(compare with ([55]) ) 



(m^ 



i^;T)) = Y. 



m^Wrr 



ie;T), 



(77) 



where 



W™(0|C;T) = (2-,5,„o)Er=o 



(n+rn p(0) n) 



(2-5,„o)E^=o 



(n+'Ti ^(0) n) 



^Er=opL(o) 



i+2E~=iE~=o 



(n+TTi p(0) n) 



i:r=oPL(o) 



(78) 



T FT ■ 

We have used the shorthands p{0\T,£_) =^ p(0) and 
W„,{I;Q\T,^) =^ W;,i(/;0) in these formulae. 



A. Pure Initial State 

We note now that, according to the first line of the 
relation ([7D|) . in the special case of a pure initial state 
the crossover 

1, if ^<^UT), 



FitT)^i _i_ I, -£^j.^ (79) 



from good to poor reversibility takes place near the crit- 
ical value ^c{T) = \j2l (rn?-)T of the strength ^ of the 
perturbation. The validity of the formula (|70p for pure 
initial states is illustrated by the numerical data plotted 
in Fig. El 

The numerical results presented in Fig. [5] imply that 
the fidelity F(^,T) decays as a function of the rever- 
sal time T the faster (approaching the exponential decay 
typical of the classical chaotic dynamics [1^) the closer 
the motion is to the semi-classical region. Within the 
interval T < the Ehrenfest time the decay is exponential 
with the rate 1/tc which describes the classical exponen- 
tial proliferation of the number of harmonics and docs 
not depend on the perturbation constant ^. 

With regard to the number of harmonics of the time- 
reversed state, expression (j78p reduces, for a pure initial 
state, (compare with Eq. (p5)) ) to 

oo 

W„M $; T) = (2 - S„,o) Y. ^n+,n{^; T) ?!;„($; T), (80) 



n=0 



where 



^„(^; T) ^ p„„(0|T, C) = |(n|e-»«"(^) |0) I' . 
In particular wo{S,\T) = F(^;T) (see Eq. (PT])). 



(81) 
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FIG. 7: (color online) Fidelity F(^; T) versus the scaled vari- 
able C/?c(r). Data correspond to: (i) h — l,go = 2; light 
pluses: T = 10, A = 0; light and black crosses: T = 50, 
A = and 25; (ii) h = 0.1, 50 = 1-5, A = 0.45 and T = 20; 
(black open circles); (iii) h — 0.01, go = 1.5, A — 0.5 and 
T = 2. (black open squares). The curves show the theoretical 
prediction of Eq. fTO)) . 



We expect again an overall exponential decay of the ex- 
citation number probability distribution: Wn>i(^;T) « 
A{£^;T) e^"^"^^"^. This assumption is well confirmed by 
numerical simulations, examples of which are presented 
in Figs. [21 O where exponential decays of both the n- 
and TTi-distributions u;„>i(^;T), >Vm(0|'^;r) is demon- 
strated in agreement with Eq. (|5D)) . It is clearly seen that 
the decay rate i'(T) does not depend on the perturbation 
strength ^. The normalization condition X)n=o '^m ^ ^ 
defines the normalization constant A{^; T) thus connect- 
ing the n-distribution with the fidelity 

Wn>i{^;T)^{l-F{i-T)) (e'^(^)-l) e-'^(^)". (82) 



A simple although a bit lengthy calculation connects, 
quite similarly to the relation expressed by Eqs. (j67p and 
([55)l . the second moment of the distribution ([50]) to the 
fluctuations of the excitation numbers: 



(m2(^;T))=2[(n2(C;T))- 

(l-F(C;T))(F(5;T)+e'- 
^ ; %■> — 






(83) 
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Notice that the ratio of the first two moments of 
the distribution ([50)1 calculated with the help of the 
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FIG. 8: (color online) Decay of Pnn(0|5,T) as a function of 
n, a.t h = 1, go = 2 and T = 50, for (a) pure (A = 0) and 
(b) mixed (A — 5) initial states. Curves from bottom to top 
correspond to different values of perturbation parameter ^ = 
^c X 2'^'^, I = —9, —6, —3, and 3. The dotted line corresponds 

to the initial distribution pnn(O) — 
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parametrization 



A(T) 



(n(e;r)) _ l-e-'^(^) 



(m2(^;r)) (fi2(e;r)) i + e--m 



(84) 



does not depend on ^ and is small if the reversal time T 
is not too small. In this case i>{T) w 2A(T) < 1. 

When the perturbation parameter ^ <C Cc(r) so that 
the Peres fidelity ((70)) is close to 1, wc obtain from Eq. 
(I83l) the ratio 



i?(e,T) 
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J( 



2 \, A(T) ; 1^ 



?e(T) 



A(Ty 
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(85) 



which characterizes the residual complexity of the re- 
versed state. This ratio remains small and therefore 
the motion remains practically reversible as long as ^ <C 
A(T). For the parameters fixed in the top panel of Fig. 8 
the ratio ^c(T)/V2A(r) sa 1.44 so that the condition of 
reversibility looks as ^ < 0.69 ^c{T) in reasonable agree- 
ment with the condition obtained above with the help of 
fidelity. 
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FIG. 9: (color online) Decay of Wm{0\£,,T) as a function of 
771, with A = and the other parameter values as in Fig. [8l 



the weighted- mean value of the lowest cumulant X2iT) 
(compare with Eq. (j67[) ) an additional positive one. Still, 
the condition ^ < (,c{T) holds as the criterion for the 
Peres fidelity to be close to one. 

Contrary to the diagonal matrix elements ([7^ , the off- 
diagonal elements of the density matrix are distorted al- 
ready in the first order approximation, 



PnS„io + i£,(n + m\h{T)\n) pn (l - -^ 



(86) 



As long as ^ is small enough only the probability of zero 
harmonic Wo is close to unity whereas the probabilities 
of other harmonics are small as ^^ : 



Wm^i(o|e;r)«e'9m(o|o,T), 

Wo{0\^;T)^l^eqo{0\0,T), 



(87) 



where 



g™^i(0|0,T) = 



^ [A+h J 



'EE^ 



T|(n-hm|n(T)|n)|^ 



dm. ^ 



|0,T) 



and 



<?o( 



|0,T)=^g,„(0|0,r). 



(89) 



B. Incoherent Initial Mixture 

When the evolution starts with a mixed initial state the 
^2-correction in the expansion ([55]) of the Peres fidelity 
contains along with the negative contribution defined by 



If the initial mixture is wide, A ^ fi,, the pre-factor 
dm ^ (l — e~3™j is small for all m < A/ft. This ex- 
plains the dip seen in Fig. [TO] (above). At the same time. 
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Fig. [To] (below) shows that the second factor fm{Q\0,T) 
is very well described by an exponential. Therefore we 
suppose that 






p,{T)m 



(90) 



To fix the two unknown functions go(0|0,T) and ii{T) 
we calculate the moments 



m('='(T) = E™=i"^'^'?™(o|o,r) = 
o;r)-^f-£)'z(A*). 



go I 



z(p) 



Then 



go(0|0;r) = -mW(T)^, 



9m(0|0;r) 



z'(/j) 



(91) 



(92) 



where the prime denotes differentiation with respect to 
the parameter fi. Besides we obtain the connection 



mW(r) _ zV) 
^ '~ m(2)(r) z"(fi) 



(93) 



which expresses the function fi{T) in the terms of the two 
lowest momenta. Under the condition A ^ 1, which is 
well justified by our numerical simulations, equation (I93p 
can be solved analytically resulting in /x(T) w 2A(T). 
Thus the two lowest moments entirely fix the probability 
distribution Wm- 

When the ratio £,/£,c{T) exceeds the unity, the dip in 
the probability distribution yV„i(0|^;T) disappears and 
the following exponential fit works well for all nonzero 
harmonics (see Fig. fTTj) : 



W™;>i- 1-Wo (e^-1) 



-fmi 



(94) 



As above we can express Wo and fi in terms of two mo- 
ments (m'^) = X^^=o"^'^^™' ^ ^ ^''^- ^ow we easily 
find that 



m 



rn'^l 



Wo = 1 - (m 

4A^ 



TO 



T^:a2 



^«l-2(m)A, 
V 4(m) A^e 



l-A 
1+A 



-2Ar. 



(95) 



where A(r) = ■^. A large number - l/2A(r) of har- 
monics have in this case similar noticeable probabilities. 
In Fig. [12] the ratio R{^,T) = {m{£„Tf)/{m^)T is 
plotted as a function of the parameter £,/^c{T). When 
this parameter is small the ratio is also small and propor- 
tional to ^2. Since the fidelity F(^; T) which describes the 
redistribution of the excitation numbers after the back- 
ward evolution is close to one under this condition, we 
conclude that the initial state is recovered with good ac- 
curacy and the motion is well reversible. On the con- 
trary, when the perturbation strength ^ exceeds the crit- 
ical value £,c{T) the fidelity becomes small and the resid- 
ual number of harmonics gets even larger than the num- 
ber of harmonics developed during the forward evolution. 







400 800 1200 

m 




1200 



m 



FIG. 10: (color online) Distribution g™ (a) and exponential 
decay of the factor fm (b), for reversal times T = 50 (bottom 
curves) and T — 100 (top curves), ai h = 1, go = 2, A = 200. 
Straight lines correspond to exp(— 2Am) [see Eq. (|53J]. 




400 800 

m 



1200 



FIG. 11: (color online) Same as in Fig. [9]but starting from a 
mixed state with A = 200. 



Therefore the evolution beconres irreversible. More pre- 
cisely, for any reversal time T, there exists an interval 
< ^ < £,c{T) of the perturbation strength ^, within 
which the quantum dynamics is approximately reversible. 
This interval, is defined by the rate of proliferation of the 
number of harmonics. In particular, this interval dinrin- 
ishes exponentially fast when the semiclassical domain is 
approached. 

At last the detailed temporal pattern is presented in 
Figs. [12] of the backward evolution for both the excitation 
number (panel (a)) as well as the number of harmonics 
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FIG. 12: Ratio R{S,, R) as a function of the parameter 
C/Cc(r), at ft = 2, ffg = 6, A = 50, and T = 50. 
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FIG. 13: (color online) Reversibility properties of quantum 
dynamics, displayed by (a) {I)t and (b) A4{t) = ^J ijn?)t , for 
different values of the perturbation parameter: from bottom 
to top, ^ = 5c X exp (//2), I = —8, . . . , —1, I = (thick black 
curve marked by the closed circle), and Z = 1, . . . , 6, at ft = 1, 
go = 2, A = 1. Circles indicate positions of the minima on 
the curves. 



VII. SUMMARY 

In this paper we have investigated the degree of stabil- 
ity and reversibility of the quantum dynamics of classi- 
cally chaotic systems beyond the semi-classical domain. 
As a measure of complexity we have used the number 
A4{t) = y/{m^)t of angular harmonics of the (initially 
isotropic) Wigner function W{I,6;t), developed during 
the evolution for the time t. This number describes the 
system's response to instantaneous rotation at that mo- 
ment by an infinitesimal angle ^ -^- 0. The number 
Ai{t) is found by calculating the distance between the 
perturbed (rotated) and unperturbed quantum distribu- 
tions. We show that, in contrast to the classical chaotic 
motion where the number of harmonics grows exponen- 
tially, A4c{t) ~ e*/"^"^ (with the rate I/tc increasing to- 
gether with the Lyapunov exponent), the number of har- 
monics of the quantum Wigner function increases, after 
the Ehrenfest time, not faster than linearly. This reveals 
much weaker sensitivity of the quantum dynamics to per- 
turbations than it is in the case of the classical dynamics. 

The relatively weak response of quantum systems to 
external perturbations makes the quantum dynamics, to 
some extent, reversible unlike the practically irreversible 
classical chaotic dynamics. To quantify this statement 
we have analyzed the degree of recovery of the initial, 
generally incoherent, mixed state after the backward evo- 
lution of the quantum distribution W{I, 9; i), rotated by 
a finite angle ^ at some reversal moment of time t = T. 
The lack of the perfect reversibility of the dynamics man- 
ifests itself by means of a redistribution of the excitation 
numbers and by the number A^(0) of harmonics of the 
Wigner function, which remains after the backward evo- 
lution. Whereas the first effect is directly described by 
the Peres fidelity the second one is, in general, revealed 
with the help of an additional infinitesimal rotation of 
the reversed state. We have shown that there exists 
a critical value £,ciT) = ^/2/M.{T) of the perturbation 
strength ^ such that the initial state is well recovered if 
^ ^ £,c{T). Reversibility disappears when the perturba- 
tion angle exceeds this value. The interval of reversibility 
< ^ < £,c{T) exponentially shrinks while approaching 
the semi-classical domain. 

Thus our analysis establishes a direct quantitative con- 
nection between the complexity of quantum phase-space 
distribution, reduced in comparison to the classical dy- 
namics, and the degree of reversibility of the quantum 
dynamics. 



(panel (b)). It is clearly seen that the time interval At 
during which the system passes in reversed order approx- 
imately the same sequence of the states, which it docs 
while evolving forward, decreases as a function of the ra- 
tio S,/£,ciT). The existence of minimal deviation (or the 
time of maximal return) during the backward evolution 
has been stressed first in UJ, UM ■ 



Ackno'wledgements 

We acknowledge support by the Cariplo Foundation 
and INFN. V.S. and O.Zh. acknowledge support by the 
RAS Joint scientific program "Nonlinear dynamics and 
Solitons" . 



15 



[1] D.L. Shepelyansky, Physica D 8, 208 (1983); G. Casati, [9] 

B.V. Chirikov, I. Guarneri, and D.L. Shepelyansky, Phys. 
Rev. Lett. 56, 2437 (1986). [10] 

[2] B.V. Chirikov, F.M. Izrailev, and D.L. Shepelyansky, [11] 
Sov. Sci. Rev. C 2, 209 (1981). 

[3] Y. Gu, Phys. Lett. A 149, 95 (1990). [12] 

[4] A.K. Pattanayak and P. Brumer, Phys. Rev. E 56, 5174 [13] 
(1997); J. Gong and P. Brumer, Phys. Rev. A 68, 062103 
(2003). [14] 

[5] K.S. Ikeda, in Quantum chaos: between order and dis- 
order, G. Casati and B.V. Chirikov (Eds.) (Cambridge 
University Press, 1995). 

[6] A. Peres, Phys. Rev. A 30, 1610 (1984). [15] 

[7] Quantity Q has been widely investigated in studies of 

the so-called quantum Loschmidt echo [g]. However, we [16] 
would like to note that, similarly to Ref. [9[ and in con- 
trast to other previous studies [3] , the backward evolution 
proceeds with the same Hamiltonian as the forward evo- [17] 
lution and the perturbation acts instantaneously only at [18] 
the reversal time T. 

[8] A review on the quantum Loschmidt echo is provided [19] 
by T. Gorin, T. Prosen, T.H. Seligman, and M. Znidaric, 
Phys. Rep. 435, 33 (2006). Ph. Jacquod and C. Petitjean, 
preprint [arXiv : 0806 ■ 0987^ -1 [quant-ph] 5 Jun 2008 



C. Petitjean, D.V. Bevilaqua, E.J. Heller, and Ph. 
Jacquod, Phys. Rev. Lett. 98, 164101 (2007). 
R.J. Glauber, Phys. Rev. 131 2766 (1963). 
G.S. Agarwal and E. Wolf, Phys. Rev. D 2, 2161 (1970); 
ibid. 2187 (1970). 

J. Schwinger, Phys. Rev. 91, 728 (1953). 
G.P. Berman, G.M. Zaslavsky, Physica A 91, 450 (1978); 
ibid. 97, 367 (1979). 

V.V. Sokolov, Nonlinear Resonance of a Quantum oscil- 
lator prep. Inst. Nucl. Phys. Siberia Div. Acad. Sci USSR 
(1978); Teor. Mat. Fiz. 61, 128 (1984) [Sov. J. Theor. 
Math. 61, 104 (1985)]. 

V.V. Sokolov, G. Benenti, and G. Casati, Phys. Rev. E 
75, 026213 (2007). 

G. Benenti and G. Casati, Phys. Rev. E 65, 066205 
(2002); G. Benenti, G. Casati, and G. Veble, Phys. Rev. 
E 67, 055202(R) (2003). 

T. Kottos and D. Cohen, Europhys Lett. 61, 431 (2003). 
Moritz Hiller, Tsampikos Kottos, Doron Cohen and Theo 
Geisel, Phys. Rev. Lett. 92 010402 (2004). 
For instance, this approach can be used for systems like 
the kicked rotor model, with the motion confined to a 
ring geometry 



